Scalar field collapse in three-dimensional AdS spacetime 
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We describe results of a numerical calculation of circularly symmetric 
■ scalar field collapse in three spacetime dimensions with negative cosmological 

. constant. The procedure uses a double null formulation of the Einstein-scalar 

, equations. We see evidence of black hole formation on first implosion of a 

scalar pulse if the initial pulse amplitude A is greater than a critical value . 
Q ' Sufficiently near criticality the apparent horizon radius rAH grows with pulse 

Q ' amplitude according to the formula rAH ~ — A^)^-^^. 



■ I. INTRODUCTION 
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^ . Dynamical gravitational collapse of matter field configurations continues to be an inter- 

esting problem in general relativity. The questions concern whether the long time static or 
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O ' stationary limit of collapse results in a black hole or a naked singularity, the details of the 

o ■ 
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u 



final outcome, and the dependence of results on matter type and symmetries of the system. 
Among the simplest such systems is the Einstein-scalar equations with minimally coupled 
Q^. massless scalar field with certain symmetries. Initial work on this system in four spacetime 
^ ! dimensions in spherical symmetry was due to Christodoulou [Q] , who showed that there exist 
classes of initial scalar configurations which give regular solutions for all time. Subsequently, 
numerical work by Goldwirth and Piran p| established that black hole formation occurs for 
^ ■ certain classes of initial data. A few years later a comprehensive study of this system by 
Choptuik answered several important questions. In particular he showed that regular 
initial data characterized by parameter p leads to a black hole with mass M {p — p^)"* 
where p^, < p is a certain critical value of p and 7 is a numerically determined exponent. 
This work also showed that the solution p = p^ is a. naked singularity. 

Since this initial work many other systems have been studied and similar behaviour found 
at the threshold of black hole formation. A semi-analytical understanding of the exponent 
in the mass formula has been obtained by assuming a critical solution and studying the 
linearized perturbation equations. For a comprehensive review and references see Ref. [Q. 

In the present work we describe a numerical simulation of minimally coupled scalar field 
collapse in three spacetime dimensions with negative cosmological constant. This system is 
of interest for a number of reasons, (i) A length scale set by the cosmological constant is 
present unlike the four dimensional Einstein-scalar case, where the only scale (apart from 
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the gravitational constant) enters via initial data, (ii) There is a static black hole solution in 
three dimensions with mass determined in units of the (negative) cosmological constant. 
In particular, there is no black hole with zero cosmological constant. (iii) The asymptotic 
spacetime is anti-deSitter (AdS) and not Minkowski, which may permit the possibility of 
reflecting an outgoing scalar field back to the interior. 

Furthermore, on a separate but related note, there is currently interest in asymptotically 
anti-deSitter metrics from string theory, where the so called AdS / CFT correspondence has 
been proposed between gravitational phenomena in the bulk of anti-deSitter spacetime and 
physical effects in certain conformally invariant Yang- Mills theories on its boundary [Q. 
The essence of the correspondence are certain proposed 'dictionary' entries which translate 
between gravitational and Yang- Mills phenomena. In particular, the correspondence relates 
low energy gravity processes with high energy Yang-Mills ones, and vice versa. Although this 
correspondence is for five dimensional AdS, there is an analog of the correspondence relevant 
to three dimensions. Thus, any critical or other behaviour found in matter collapse in the 
bulk AdS spacetime may provide hints for a dictionary entry for the AdS/CFT conjecture. 
As such results of the present work are potentially of interest for string theory. 

There has been other work, both analytical and numerical on matter collapse in three 
dimensions. There are Vaidya-like solutions |0] for collapsing dust, with and without pres- 
sure. Analytic results have been found for the collapse of dust shells IQ and point particles 
||lO| , p]T|| . A recent numerical study of this same system has been carried out by Pretorius 
and Choptuik using a different coordinatization and numerical technique . More recently 



an exact critical self-similar solution has been found [0, and the nature of the singularity 



within the black hole has been studied using a model system [|1^ . 

This paper is organized as follows. The next section presents the system of equations and 
the details of the boundary conditions we use. Section III describes the numerical procedure. 
The final section describes the results and outline of possible future work. 



II. EINSTEIN-SCALAR EQUATIONS 

The Einstein-scalar equations with cosmological constant A = —1//^ are 
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Gab - ^ gab = 2vrTafe, (l) 



Tab = da(j)db(j) - ^gab9^'^dc(j)dd(f). (2) 
These may be rewritten in the form 

2 

Rab + gab = 27rda4>db(j). (3) 

The general circularly symmetric metric in three spacetime dimensions in double null 
coordinates {u, v) may be written as 

ds"^ = dab dx"'dx^ = —21 g{u, v) r'{u, v) dudv + r'^{u, v) dO"^, (4) 
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where we are taking g and r to be dimensionless functions, {u, v) and / have dimensions of 
length, and /' = df/dv, f = df/du. In this parametrization, coordinate freedom has been 
fixed since there are only two metric functions to be determined. Time evolution may be 
viewed as occur ing with respect to the coordinate u. 

The vv component of the Einstein equations is the constraint 

— = 2v^(</>')^ (5) 
gr 

on a fixed u surface. The 66 component is the time (or u) evolution of r 

Ir' = —rr'g. (6) 

Finally, the time evolution of the scalar field is provided by the Klein-Gordon equation 

(r0') + (r0)' = 0. (7) 

The remaining components of the Einstein equation are identically satisfied as consequence 
of these three equations. 

If = these equations admit a black hole solution first discovered in . In the standard 
(r, t) coordinates the black hole metric is 



-1 



ds' = -\^[j^ -Mjdi'+y[j^ -Mj rfr^ + rW. (8) 

AdS space is recovered by setting M = — 1. We can label the outgoing radial null geodesies 
of the AdS metric by u and the ingoing ones by v. These coordinates are defined by 

u = t~l tan"^ (^0 V = t + l tan"^ (9) 

Therefore AdS space in double null coordinates with our choice of dimensions and 
parametrization is given by 

r{u,v) = tan , 9iu,v) = 1. (10) 

Similarly, the BTZ black hole in the corresponding null coordinates is given by Q 



r{u, v) = y M tanh 



21 



[u — v] 



g{u,v) = l. (11) 



Thus r = corresponds to u = v. In the following we will use this as one of the boundary 
conditions for the full set of coupled equations. 

Equations (§-|^) may be rewritten in a first order initial value form suitable for numerical 
integration as follows: On a constant u surface, integrating the constraint equation (H) gives 



^Note that this is the BTZ solution for r < \fM. For r > \/M the tanh is replaced by coth. 
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(12) 



g = ki{u) exp 

where ki{u) is an arbitrary function of u. 
Integrating (D gives 

r = -y£rf5 (^rr') + ^, (13) 

where k2{u) is an arbitrary function of u. Note that in both of these equations we are 
integrating from m to f which is to correspond to integrating outward from r = 0. 
Finally, using r from (|T3|), and the definition 

the Klein-Gordon equation (0) gives 
where 

g(u,v) = / dv {grr') . (16) 

J u 

The definition of h Eqn. ([I^ ) gives 



The function /c3('u) is fixed by noting from Eqn. ( |T4|) if = C then h = C, where C is a 
constant. Comparing this with Eqn. ( |T7| ) gives /i;3(M) = 0. 

The equations (|12D, ([T3|) , and ([T5| ) give the required null initial value formulation for the 
three-dimensional Einstein-scalar system with negative cosmo logical constant A = — 1//^. 
These equations contain two arbitrary functions ki{u) and k2{u), which are to be fixed by 
boundary conditions. With these functions fixed, initial data sets r{u = 0, v) and h{u = 0, v) 
(derived from 0) can be evolved to later times u. 

Let us now see how the arbitrary functions in the equations are fixed. The most natural 
boundary condition is the requirement that the metric is regular at the origin r = 0, which 
is a timelike line. In our formulation r{u, v) is a field variable. Therefore we must first define 
where r = is in order to impose regularity conditions there. As noted above the natural 
choice motivated by AdS and the BTZ black hole is r = at u = f , i.e. 

r{u,u) = 0, (18) 

which implies 

f{u,u) = —r'{u,u). (19) 
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With this choice for r = 0, there are two possible boundary conditions that may be imposed 
at r = 0. Before describing these we set ki{u) = 1 for simphcity. 
(i) No conical singularity at r = 0. This gives the condition 



r{u,u) = —n^g{u,u)/2l. 
where n is an integer. (See the Appendix.) Using this in Eqn. ( |T3| ) gives 



koiu) 



2 ~ "Y' 



(20) 



(21) 



which is to be used in the evolution equations (|r^) and (|T^. 
(ii) The uv metric function is unity: 



21 g{u,u) r'{u,u) = —2lg{u,u) r{u,u) = 1. 



Therefore 



r[u, u 



2lg{u,u)' 



and eqn. (O) gives 



1 
2 



(22) 



(23) 



(24) 



We shall use k2{u) = — 1/2. In summary the first order equations to be integrated numeri- 
cally are 



^ I 21 



where 



g = exp 
g = dv (grr' 



f/:-(7<" 



1 r ,^ f h r'' 

dv 



2y^ Ju \Vr ) ' 

Notice that with h = (p = our equations with k2 an arbitrary constant reduce to 



r M, V 



r2 - 2k2 
21 



which is to be integrated with the initial condition r{u,u) = 0. The solution is 



riu, v) = \/2k2 tanh 



21 



'-{u — v) 



(25) 
(26) 

(27) 
(28) 
(29) 

(30) 
(31) 
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Comparing with (|Tl|), we see that 

M = 2k2. 



(32) 



Thus the value k2 = —1/2, which we use for numerical integration, corresponds to M = — 1, 
which is 3- dimensional AdS. Notice that although this is the natural case to study in that 
it is most similar to the spherically symmetric collapse in four dimensions (because r = 
is a regular point initially), the other choices for k2 corresponding to M = 0, M > and 
— 1 < M < are also interesting. 



III. NUMERICAL METHOD 

Our numerical scheme is based on the approach used in and its refinements in [0]. 
The evolution equations (p5|-p6D are converted into set of coupled ordinary differential 
equations on a grid by the correspondence 

h{u,v) ^ hn{u), r{u,v) ^ rn{u), (33) 

and similarly for the derived variables g{u,v), g{u,v), and (p{u,v) where n = 0,---,N 
specifies the v grid. The differential equations are therefore 



9n _ ]_ 
' I 21 



(34) 



hn = J {K - 0n) \9n~^^- ^ j (35) 

We take the initial value functions to be 

r(0,w)=w, (36) 
with /i(0, v) determined by this and the initial scalar field configuration 

<i){u = 0, r) = Ar^ exp (-^^^) , (37) 

where A, a and rg are initial data parameters of the scalar field. The initial specification of 
the other functions is completed by computing the integrals for gn and gn using Simpson's 
rule. 

The boundary conditions at fixed u are 

Tk = 0, gk = 0, gk = 1. (38) 

where k is the index corresponding to the position of the origin r = 0. (In the algorithm 
used all grid points < i < k — 1 correspond to ingoing rays that have reached the origin 
and are dropped from the grid; see below). These conditions are equivalent to r{u,u) = 0, 
g\r=o = g{u,u) = 1 etc., and guarantee regularity of the metric at r = 0. Notice that for 
our initial data, (pk and hence hk are initially zero, and therefore remain zero at the origin 
because of Eqn. (|35|) . 
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Our numerical procedure is similar to that used in Refs. . Time evolution with respect 
to the coordinate u is initiated via a fourth order Runge-Kutta scheme. The constraint 
integrals are recalculated on each u slice. As evolution proceeds the entries in the r„ array 
sequentially reach r = 0, at which point they are dropped from the grid. Thus the radial 
grid loses points with evolution. This is similar to the procedure used in and ||^. At 
each u step, a check is made to see if an apparent horizon has formed. This is done in the 
standard way by calculating and plotting the function 

g-%rd,r = (39) 

Evolution terminates if either an apparent horizon forms, or the scalar field pulse reflects 
off r = 0. In practice we stop the evolution if \dur\ becomes 10~^. 

To improve numerical accuracy near r = we follow a procedure similar to [0], where all 
functions on a constant u surface are expanded in power series in r at r = 0, and the first 
three values of the constraint integrals are calculated using the respective power series. We 
write 



h = ho + hir (40) 
from which the following relations are derived: 

hir TiK^r"^ r^lngf , , 

<P = h, + ^, ^ = l + _^, ^ = _ + _y. (41) 

In practice, because of our boundary conditions, turns out to be of order 10~^. The 
remaining — 3 values of these functions are computed using Simpson's rule for equally 
spaced points. This substantially improves accuracy near r = 0. f— derivatives of functions 
are calculated using // = (/j+i — fi-i)/2Av with end point values determined by linear 
extrapolation: f[ = 2/^ - f!^ and f'j^ = 2f'j^^^ - f'^_^. 

All numerical calculations were performed for / = 1 with up to 1400 v-grid points and 
scalar field parameters a = 0.3 and ro = 1.0. The code is written in C and was executed on 
a Sun workstation. 

Before describing our results in the next section, we comment further on the numerical 
procedure outlined above. 

Firstly, because we use double null coordinates with the number of v grid points de- 
creasing as ingoing null geodesies cross r = 0, we are not able to follow a reflected pulse 
back out toward infinity. This is likely a shortcoming because after first implosion some 
percentage of the scalar field intensity will be reflected back toward infinity and reach it 
in finite time (due to the compactness of AdS space). It is then possible that some part 
of the field would be reflected back toward the origin from the boundary at infinity. This 
would tend to increase the mass of the black hole. Multiple reflections could occur several 
times until a static limit is reached. Our procedure is limited to evaluating the mass of the 
black hole on first implosion of the scalar field pulse. We can assume however that the mass 
increase due to multiple reflections would be small compared to the mass obtained on first 
implosion. 
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Secondly, again due to the loss of points from the grid, we are not able to observe the 
behaviour of the scalar field very near criticality. However, the calculations performed in 
indicate that there is continuous self similarity in the near critical regime. f\ 



IV. RESULTS AND COMMENTS 

After several simulations for different values of the initial amplitude A of the scalar field, 
we arrive at the following conclusions. For amplitudes smaller than A^.. ~ 0.1248 the field is 
reflected back in totality from the origin and no black hole is formed. For A > A^ however, 
we observe black hole formation. 
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FIG. 1. Logarithmic plot of the apparent horizon radius vah versus the initial amplitude of the 
imploding scalar field (A — A^,). 

Figure 1 shows a logarithmic plot of the apparent horizon radius r^H versus the initial 
amplitude A. The plot shows that for A sufficiently close to A^ the behaviour is nearly 
linear (the first six points). This is the region of interest for us since it appears to exhibit 
the type of scaling behaviour observed in many other collapse studies We ignore the 

other part of the plot which shows some oscillatory behaviour and eventually a flattening of 



^ As in our case, and for similar reasons, this work also focusses on a search for critical and scaling 
behaviour on first implosion of the scalar pulse. 
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the curve, because this is further from criticahty. We thus estimate the value of the critical 
exponent by considering the linear part. This yields the result 



TAH -{A- A^y (42) 

where 7 ~ 0.81. Given that black-hole mass is proportional to for the BTZ black- 
hole, our mass exponent estimate is 1.62. Pretorius and Choptuik [|T2| used a different 
method proposed by Garfinkle to estimate the critical exponent. They studied the subcritical 
[A < A^) scaling behaviour of the maximum value of the scalar curvature R to obtain a 
value for 7 in the range 1.15 — 1.25. 

There are several directions for future study using the double null procedure outlined 
here. These include computations for other values of the constant k2 and other initial scalar 
configurations. In addition, the self similar behaviour of the scalar field near criticality found 
numerically in and analytically in |TB[, can likely be reproduced using the method of 



grid point restoration used in 0. 

In connection with the AdS/CFT conjecture, it will be useful to carry out similar cal- 
culation for 5-dimensional AdS space. The computational procedure used here requires 
relatively minor modifications to study that case. 
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V. APPENDIX 

The condition of no conical singularity at the origin r = is derived. The parallel 
transport of a vector A" on a closed loop around the origin should be such that this vector 
comes back to itself after a complete turn. The condition for parallel transport of A°- along 
a curve tangential to is 

eVaA^ = 0. (43) 

Moving on a curve of constant u and v, only 6 varies from to 27r, i.e. t"" = {d/dOY and 
the above equation becomes 

VeA^ = deA^ + T'^g^A' = 0. (44) 

The relevant non-zero F are 

^ee = - Tl, = - (45) 

rL = K. = (46) 
The parallel transport equation is then 
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(47) 



Diagonalizing the matrix on the right we get 3 eigenvalues A = 0, ~- Iii the eigenvectors 
basis we would thus have the following differential equation to solve for each component 

(9eA" = A(°M". (48) 

This is solved by A°'(6) = y4"(0) exp (^A'-"^^^ . Requiring that A" be unchanged after a 
complete non-trivial turn around the origin requires 



A"(27r) = A"(0) v4"(0)exp (27rA(")) = A°(0) A^^^ 



(49) 



where n eZ. The simplest non-trivial case is A*^"-* = z, and this becomes a boundary condition 
if we look at the case where a is the third eigenvalue : 



at 



0. 



(50) 



This first boundary condition is thus 
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' = -2 



at 



r = 0. 



(51) 
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